Modeling the response of Norway spruce tree-ring carbon and oxygen isotopes to selection harvest on a drained peatland forest.

Continuous cover forestry (CCF) has gained interest as an alternative to even-aged management particularly on drained peatland forests. However, relatively little is known about the physiological response of suppressed trees when larger trees are removed as a part of CCF practices. Consequently, studies concentrating on process-level modeling of the response of trees to selection harvesting are also rare. Here, we compared, modeled and measured harvest response of previously suppressed Norway spruce (Picea abies) trees to a selection harvest. We quantified the harvest response by collecting Norway spruce tree-ring samples in a drained peatland forest site and measuring the change in stable carbon and oxygen isotopic ratios of wood formed during 2010-20, including five post-harvest years. The measured isotopic ratios were compared with ecosystem-level process model predictions for ${\kern0em }^{13}$C discrimination and ${\kern0em }^{18}$O leaf water enrichment. We found that the model predicted similar but lower harvest response than the measurements. Furthermore, accounting for mesophyll conductance was important for capturing the variation in ${\kern0em }^{13}$C discrimination. In addition, we performed sensitivity analysis on the model, which suggests that the modeled ${\kern0em }^{13}$C discrimination is sensitive to parameters related to CO2 transport through stomata to the mesophyll.

Potential of continuous cover forestry on drained peatlands to increase the carbon sink in Finland.

Land-based mitigation measures are needed to achieve climate targets. One option is the mitigation of currently high greenhouse gas (GHG) emissions of nutrient-rich drained peatland forest soils. Continuous cover forestry (CCF) has been proposed as a measure to manage this GHG emission source; however, its emission reduction potential and impact on timber production at regional and national scales have not been quantified. To quantify the potential emission reduction, we simulated four management scenarios for Finnish forests: (i) The replacement of clear-cutting by selection harvesting on nutrient-rich drained peatlands (CCF) and (ii) the current forest management regime (BAU), and both at two harvest levels, namely (i) the mean annual harvesting (2016-2018) and (ii) the maximum sustainable yield. The simulations were conducted at the stand scale with a forest simulator (MELA) coupled with a hydrological model (SpaFHy), soil C model (Yasso07) and empirical GHG exchange models. Simulations showed that the management scenario that avoided clear-cutting on nutrient-rich drained peatlands (i.e. CCF) produced approximately 1 Tg CO2 eq. higher carbon sinks annually compared with BAU at equal harvest level for Finland. This emission reduction can be attributed to the maintenance of a higher biomass sink and to the mitigation of soil emissions from nutrient-rich drained peatland sites.

Soil GHG dynamics after water level rise - Impacts of selection harvesting in peatland forests.

Managed boreal peatlands are widespread and economically important, but they are a large source of greenhouse gases (GHGs). Peatland GHG emissions are related to soil water-table level (WT), which controls the vertical distribution of aerobic and anaerobic processes and, consequently, sinks and sources of GHGs in soils. On forested peatlands, selection harvesting reduces stand evapotranspiration and it has been suggested that the resulting WT rise decreases soil net emissions, while the tree growth is maintained. We monitored soil concentrations of CO2, CH4, N2O and O2 by depth down to 80 cm, and CO2 and CH4 fluxes from soil in two nutrient-rich Norway spruce dominated peatlands in Southern Finland to examine the responses of soil GHG dynamics to WT rise. Selection harvesting raised WT by 14 cm on both sites, on average, mean WTs of the monitoring period being 73 cm for unharvested control and 59 cm for selection harvest. All soil gas concentrations were associated with proximity to WT. Both CH4 and CO2 showed remarkable vertical concentration gradients, with high values in the deepest layer, likely due to slow gas transfer in wet peat. CH4 was efficiently consumed in peat layers near and above WT where it reached sub-atmospheric concentrations, indicating sustained oxidation of CH4 from both atmospheric and deeper soil origins also after harvesting. Based on soil gas concentration data, surface peat (top 25/30 cm layer) contributed most to the soil-atmosphere CO2 fluxes and harvesting slightly increased the CO2 source in deeper soil (below 45/50 cm), which could explain the small CO2 flux differences between treatments. N2O production occurred above WT, and it was unaffected by harvesting. Overall, the WT rise obtained with selection harvesting was not sufficient to reduce soil GHG emissions, but additional hydrological regulation would have been needed.

Estimating intra-seasonal photosynthetic discrimination and water use efficiency using δ13C of leaf sucrose in Scots pine.

Sucrose has a unique role in recording environmental and physiological signals during photosynthesis in its carbon isotope composition (δ13C) and transport of the signal to tree rings. Yet, instead of sucrose, total organic matter (TOM) or water-soluble carbohydrates (WSC) are typically analysed in studies that follow δ13C signals within trees. To study how the choice of organic material may bias the interpretation of δ13C records, we used mature field-grown Scots pine (Pinus sylvestris) to compare for the first time δ13C of different leaf carbon pools with δ13C of assimilates estimated by a chamber-Picarro system (δ13CA_Picarro), and a photosynthetic discrimination model (δ13CA_model). Compared with sucrose, the other tested carbon pools, such as TOM and WSC, poorly recorded the seasonal trends or absolute values of δ13CA_Picarro and δ13CA_model. Consequently, in comparison with the other carbon pools, sucrose δ13C was superior for reconstructing changes in intrinsic water use efficiency (iWUE), agreeing in both absolute values and intra-seasonal variations with iWUE estimated from gas exchange. Thus, deriving iWUE and environmental signals from δ13C of bulk organic matter can lead to misinterpretation. Our findings underscore the advantage of using sucrose δ13C to understand plant physiological responses in depth.

Explicitly accounting for needle sugar pool size crucial for predicting intra-seasonal dynamics of needle carbohydrates δ18 O and δ13 C.

We explore needle sugar isotopic compositions (δ18 O and δ13 C) in boreal Scots pine (Pinus sylvestris) over two growing seasons. A leaf-level dynamic model driven by environmental conditions and based on current understanding of isotope fractionation processes was built to predict δ18 O and δ13 C of two hierarchical needle carbohydrate pools, accounting for the needle sugar pool size and the presence of an invariant pinitol pool. Model results agreed well with observed needle water δ18 O, δ18 O and δ13 C of needle water-soluble carbohydrates (sugars + pinitol), and needle sugar δ13 C (R2 = 0.95, 0.84, 0.60, 0.73, respectively). Relative humidity (RH) and intercellular to ambient CO2 concentration ratio (Ci /Ca ) were the dominant drivers of δ18 O and δ13 C variability, respectively. However, the variability of needle sugar δ18 O and δ13 C was reduced on diel and intra-seasonal timescales, compared to predictions based on instantaneous RH and Ci /Ca , due to the large needle sugar pool, which caused the signal formation period to vary seasonally from 2 d to more than 5 d. Furthermore, accounting for a temperature-sensitive biochemical 18 O-fractionation factor and mesophyll resistance in 13 C-discrimination were critical. Interpreting leaf-level isotopic signals requires understanding on time integration caused by mixing in the needle sugar pool.

Does growing atmospheric CO2 explain increasing carbon sink in a boreal coniferous forest?

The terrestrial net ecosystem productivity (NEP) has increased during the past three decades, but the mechanisms responsible are still unclear. We analyzed 17 years (2001-2017) of eddy-covariance measurements of NEP, evapotranspiration (ET) and light and water use efficiency from a boreal coniferous forest in Southern Finland for trends and inter-annual variability (IAV). The forest was a mean annual carbon sink (252 [ ± 42] gC m-2a-1 ), and NEP increased at rate +6.4-7.0 gC m-2a-1 (or ca. +2.5% a-1 ) during the period. This was attributed to the increasing gross-primary productivity GPP and occurred without detectable change in ET. The start of annual carbon uptake period was advanced by 0.7 d a-1 , and increase in GPP and NEP outside the main growing season contributed ca. one-third and one-fourth of the annual trend, respectively. Meteorological factors were responsible for the IAV of fluxes but did not explain the long-term trends. The growing season GPP trend was strongest in ample light during the peak growing season. Using a multi-layer ecosystem model, we showed that direct CO2 fertilization effect diminishes when moving from leaf to ecosystem, and only 30-40% of the observed ecosystem GPP increase could be attributed to CO2 . The increasing trend in leaf-area index (LAI), stimulated by forest thinning in 2002, was the main driver of the enhanced GPP and NEP of the mid-rotation managed forest. It also compensated for the decrease of mean leaf stomatal conductance with increasing CO2 and LAI, explaining the apparent proportionality between observed GPP and CO2 trends. The results emphasize that attributing trends to their physical and physiological drivers is challenged by strong IAV, and uncertainty of LAI and species composition changes due to the dynamic flux footprint. The results enlighten the underlying mechanisms responsible for the increasing terrestrial carbon uptake in the boreal zone.